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Abstract. The azimuthal substructure observed in some debris disks, as exemplified by e Eridani, is usually attributed to 
resonances with embedded planets. In a standard scenario, the Poynting-Robertson force, possibly enhanced by the stellar wind 
drag, is responsible for the delivery of dust from outer regions of the disk to locations of external mean-motion planetary 
resonances; the captured particles then create characteristic "clumps". Alternatively, it has been suggested that the observed 
features in systems like £ Eri may stem from populations of planetesimals that have been captured in resonances with the 
planet, such as Plutinos and Trojans in the solar system. A large fraction of dust produced by these bodies would stay locked 
in the same resonance, creating the dusty clumps. To investigate both scenarios and their applicability limits for a wide range 
of stars, planets, disk densities, and planetesimal families we construct simple analytic models for both scenarios. In particular, 
we show that the first scenario works for disks with the pole-on optical depths below about ~ 10~ 4 -10~ 5 . Above this optical 
depth level, the first scenario will generate a narrow resonant ring with a hardly visible azimuthal structure, rather than clumps. 
It is slightly more efficient for more luminous/massive stars, more massive planets, and planets with smaller orbital radii, but 
all these dependencies are weak. The efficiency of the second scenario is proportional to the mass of the resonant planetesimal 
family, as example, a family with a total mass of ~ 0.01 to 0.1 Earth masses could be sufficient to account for the clumps of 
e Eridani. The brightness of the clumps produced by the second scenario increases with the decreasing luminosity of the star, 
increasing planetary mass, and decreasing orbital radius of the planet. All these dependencies are much stronger than in the 
first scenario. Models of the second scenario are quantitatively more uncertain than those of the first one, because they are very 
sensitive to poorly known properties of the collisional grinding process. 
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1 . Statement of the problem 

Various kinds of structures — offsets, wing asymmetries, 
warps, voids, clumps, rings, spirals, inner gaps — are ei- 
ther seen in the images of resolved debris disks or have 
been retrieved from analyses of spectral energy distributions. 
Accordingly, a variety of mechanisms h ave been sugges ted to 
explain the structure formation (see, e.g. Augereau 1 2004, for a 
recent review). Embedded planets can create stru ctures through 
secular (long-term) perturbations ( Mouillet et al. 1997), by res- 

as well as by gravita- 



secular (long-term) perturbations ( Mouil 
onant perturbations (Liou & Zook 1999) 



tional scattering of the disk particles ( Moro-Martin & Malhotra 
200s |). Some featu res may stem from recent stellar flybys 
llKalas et alJ l200ll) : others may reflect recent major colli- 
sional events between large planetesimals (fWvatt & Dentl2002t 
iKenvon & Bromle\l2005t iGrigorieva et all20 06 ) . 

Here we address the azimuthal substructure seen in sev- 
eral disks, best exemplified by the 0.85 Gyr-old K2V star, 
e Eridani. Submillimeter images obtained with JCMT/SCUBA 
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sighted a few bright blobs or clumps in this nearly pole- 
on disk (Greaves et al. 1998). New observations of the same 
disk jGreaves et alJl2005l) revealed signs of the disk rotation, 
demonstrating that three out of six clumps are real and indi- 
cating that the observed rotation rate is compatible with the 
Keplerian velocity at the location of the clumps. A similar 
clumpy structure has recently been found i n the pole-on disk 
around another 1 Gyr-old K star, HD 53143 jKalas et alJ2006l) . 
Clumps have also been reported for other disks seen at moder- 
ate inclin ations, most notably Fomalhaut Jljolland et aljE()03l) 
and Vega dHolland et al.ll998HWilner et all2002l) . 

It is natural to attribute the formation of blobs to the 
resonant perturbations exerted on dust grains by a planet. 
But how can the resonances cr eate the clum ps? One mech- 
anism, originally proposed by Gold] i 19751) and discussed 
very intensively in the past decade, which we call sce- 
nario I, can be summarized as follows. Dust particles in di- 
lute debris disks are known to steadily migrate inwards due 
to the Poynting-Robertson (P- R) force, possibly enhanced 
by the stellar wind (e.g. [W eidenschilling & Jackson] Il993t 
iBeauge & Ferraz-Melldl 19941 Liou & Zooklll997l) .lfthere is 
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a planet in the disk, the particles reach locations of exte- 
rior mean-motion resonances (MMRs) with the planet and get 
trapped, yielding characteristic density patterns. 



iLiou et alJ ( fl99^) and lLiou & ZookH 19991) applied this idea 
to the solar system. By modeling perturbations induced by 
jovian planets on the Edgeworth-Kuiper Belt dust disk, they 
found efficient resonant trapping of dust by Neptune (which 
produces arcs of dust co-orbital with the planet) and efficient 
ejection of dust out of the solar system by Jupiter and Saturn. 
Were the solar system observed from outside, the presence of 
at least Neptune and Jupiter would be obvious merely from an- 
alyzing images of the dust disk. While the large-scale clumps 
on the outskirts of the solar system still escape observational 
detection, there is one, so far the only, case where scenario I 
is observed at work: the asymmetric reson ant ring of aster- 
oidal dust around the Earth orbit. Predicted by Jackson & Zook 
il989fc . it was identi fied in IRAS (foermott et al] 1 19941) and 
COBE/DIRBE data jReach et alJ Il995l) . Search of a simi- 
lar ring aroun d the Mars orbit was not successful, however 
jKuchner et al.l200C4) . 



In the case of e Eri. ILiou et all J2000h- lOzernov etal. 

1 2000I) . lOuillen & Thorndikd J2002I) . and beller & Maddisoii 
(2005J) have tried to find particular orbital parameters and 
masses of planet, or planets, that may reproduce the observed 
substructure in the disk. A systematic overview of resonant 
structures that a single planet can form in a debris disk through 
this mechanism was painted by Kuchner & Holman (2003). 



A major problem with this mechanism is that catas- 
trophic grain-grain collisions may smear out the planet- 
induced structure already at moderate optical depths, because 
the collisional lifetimes appear to be shorter than the P- 
R timescal e that determines the dust injecti o n rate into the 
reson a nce JLecavelier des Etangs et all 1 19961: lLagrange et alJ 
2000). IWvattl d2003l " 2006) suggest an alternative scenario, 
hereinafter scenario II, in which the clumps stem from plan- 
etesimals that have been captured in planetary resonances, pos- 
sibly as a result of the planet's migration. Such resonant pop- 
ulations have long been known in the solar system — Plutinos 
are locked in the 3:2 MMR with Neptune, while Greeks and 
Trojans reside in the primary resonance with Jupiter. A large 
fraction of the dust which these bodies produce may stay 
locked in the same resonance and could well account for the 
observed clumps. 

In this paper, we make a comparative analysis of both sce- 
narios and estimate typical timescales and efficiency of both 
scenarios with simple analytic models. The goal is to investi- 
gate and compare the efficiency of both scenarios for a wide 
range of essential parameters, such as the mass of the central 
star, mass and orbital radius of the embedded planet, as well 
as the optical depth of the "background" disk (in the first sce- 
nario) and properties of the resonant family of planetesimals (in 
the second one). Scenarios I and II are considered in Sects. 2 
and 3, respectively. Section 4 contains our conclusions and a 
discussion. 



2. Scenario I 

We start with an analysis of scenario I, in which dust particles 
are brought from afar to resonance locations by the P-R effect 
and are captured in resonances. A resonance implies a (p + 
q) : p-commensurability of the mean motions of the planet and 
the grain. If a p is the semimajor axis of the planet orbit, the 
resonant value of the grain's orbital semimajor axis is 



!fl ,(l-P)W£±i 



2/3 



(1) 



with p being the radiation pressure to gravity ratio for the dust 
grain (assumed to be a perfect absorber): 



where L*/L© and M*/Mq are luminosity and mass of the star 
in solar units and s is the grain radius. In what follows, we 
treat as a parameter and calculate luminosity by using the 
standard L* °< Aff relation for main-sequence stars. Throughout 
the paper, the bulk density of 2gcm~ 3 is adopted. To keep 
our treatment simple, we confine our analysis to a circular 
planetary orbit. The typical resonan t dynamics can then be 
summarized as follows (e.g.lWeiden schillin^ & Jacks on! 1993t 
iBeauge & Ferraz-Mellol Il994l ILiou & Zookl Il997l) . When 
trapped, a dust grain preserves the semimajor axis. If A, and 
Xp are longitudes of the grain and planet, respectively, and 6) 
is the longitude of pericenter of the grain orbit, the resonant 
argument 



<E> = (p + q)X - pX P - q&> 



(3) 



libr ates with an amplitude A around the v alue w 180° (cf. 
iKuchner & Holmanl2003llWvatl2003ll2006l) . T he eccentricity 
gradually pumps up. For initially circular orbits JLiou & Zookl 
Il997l) 
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where 



B 



2GM*|3 



(4) 



(5) 



GM* standing for the gravitational parameter of the central star 
and c for the speed of light. When the eccentricity reaches a 
certain value e res , the particle becomes a planet-crosser and 
shortly thereafter either collides with the planet or gets ejected 
out of the resonance zone after a close encounter. The orbital 
inclination ; exponentially decreases on the same timescale as 
e increases; here we assume < i < £, where e is a (small) 
semi-opening angle of the disk. 

Note that some grains that used to be trapped in a reso- 
nance may then continue migrating toward the star, instead of 
being ejected from the system. This is, however, more typical 
of resonances with very low-mass planets — e.g. in the Earth 
resonant ring mentioned above. For such planets, overlapping 
of first-order resonances and the onset of chaos, as well as ex- 
ternal perturbations by other planets, may also play a role, mak - 
ing the dynamics more complicated ( Marzari & Vanzanill9 94). 
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Since dust structures induced by terrestrial planets around other 
stars are far beyond observational limits, our analysis is con- 
fined to at least Neptune/Uranus-mass planets and thus ignores 
any temporary resonant capture of particles. 

To construct a simple model, we choose a "typical" grain 
radius s and a corresponding p-value. This is the size of the 
grains that compromises the cross section-dominating size and 
the size for which the probability of capture in the resonance 
{pres below) becomes significant. The former corresponds to 
3 of several tenth s jKrivov et alJl200Ct iThebault et al]l2003t 
Krivov et al. 2006). The latter requires p < 0. 1 (e.g. Liou et al. 
19961 Liou & Zookl 19991) : the smaller p\ the larger p res . In our 
model, we treat 3 as a free parameter and use P = 0. 1 as a de- 
fault value in numerical examples. 

For a certain MMR and a given P ratio, the res- 
onant value of semimajor axis a is given by Eq. Q. 
The particles drift towards the star by P-R drag. Assume 
that h + grains cross r — a per unit time. Let p res be 
the probability of trapping in the resonance. Its value de- 
pends on size, initial eccentricity, and the inclination of 

the grain orbit, planet's mass M n , and varies from one 

I — ~ II 1 

resonance to another (e.g. Weidenschilling & Jackson 1993; 

Beauge & Ferraz-Mellolll994 iLiou et alJll996t iLiou & Z ook 
19971 For instance. Liou & Zookl dl999h report that p res 
increases from about 0.2 for P = 0.4 to about 0.5 for 
P = 0.05 for the Edgeworth-Kuiper belt particles' trap- 
ping by Neptune. Inside a certain range of these pa- 
rameters, however, p res is not very sensitive to their 
variation. We performed a number of numerical integra- 
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is very well known (e. 


? 

g. Weidenschilling & Jackson 


1993 


Beauge & Ferraz-Mello 


1 19941 lLiou& Zookl 119971 


Ouillen & Thorndike 2002; Deller & Maddison 2005). For the 



parameters of the e Eri system, with initial eccentricities ^ 0. 1, 
inclinations below 10°, and a planet at 40 AU with 0.1 Jupiter 
mass, the trapping probability is between 10% and almost 
100% for many resonances (2:1, 3:2, 8:5, 5:3, 9:5, and oth- 
ers). In our model, p res is a parameter. In numerical examples, 
P res = 0.5 is used as a default value. 

Besides p res , another quantity we need is e res , the value 
of orbital eccentricity that the captured grains have to develop 
to quit the resonance. On the basis of the papers cited above 
and our own numerical integrations, the typical values lie be- 
tween 0.1 and 0.4. Like P and p res , e res is a model param- 
eter. As a "standard" value, we adopt e res = 0.2. Note that 
grains do not necessarily reach e = e res . Let us denote by T res 
the time interval in which e res is attained. As the particles are 
lost not only because they pump up the orbital eccentricities 
to the resonance-quitting value e res , but also are eliminated by 
mutual collisions, their mean actual lifetime T is shorter than 
T res . Therefore, their actual maximum eccentricity e max , i.e. the 
maximum eccentricity achieved by the particles in their life- 
time T, will be smaller than e res . 

We now introduce an annulus in the disk around the loca- 
tion of a given MMR (Fig.G), defined as follows. We require 
that the annulus width matches the maximum radial excursion 
of the trapped particles: from a{\ +e max ) to a (I — e max ). At any 
instant in time, the annulus contains two populations of parti- 



Clump 




Grain orbit 



Annulus 



Fig. 1. Schematic of the resonant and non-resonant compo- 
nents of the dust disk in the vicinity of a given resonance. 
The particular resonant population with two clumps shown here 
corresponds to the 3:2-MMR with emax = 0.2andA = 20°. 

cles. One is a non-resonant "background" population of parti- 
cles, including those that have not reached the resonance loca- 
tion yet, as well as those that drifted through it without getting 
trapped. It might also include grains that continue drifting to- 
ward the star after quitting resonance. If such grains are present 
in sufficient amounts, p res should be re-defined to account for 
them. All these particles have a uniform azimuthal distribution. 
Another population is an unevenly distributed (clumpy) popu- 
lation of resonant grains. Hereinafter, all quantities with sub- 
script "0" characterize the background population and those 
without subscript the resonant population. For instance, iiq and 
n stand for the total number of background grains and resonant 
grains in the annulus, respectively. 

We now construct a simple kinetic model. By applying any 
standard method of the kinetic theory, or simply by "counting" 
the particles injected into and destroyed in the annulus per unit 
time, one finds the balance equation: 



dn 
dt 



— = n' p r 



n 

Trf>\ 



'coll 



n 

Tcoll ' 



(6) 



where Ts are timescales as explained below. As n does not have 
a grain size/mass as argument, Eq. can be referred to as the 
Boltzmann equation. Since n does not have any grain velocity 
argument eith er, it can also be term ed the Smoluchowski equa- 
tion. See, e.g.. lKrivov et al.N2005l) and references therein for a 
detailed discussion of a general form of the balance equation. 

The resonant population implies several natural timescales 
associated with the terms on the right-hand side of Eq. ©. The 
first is Tpr, the P-R time to drift from a ( 1 + e max ) to a ( 1 — emax) . 
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The second is T res , the time of eccentricity pumping to the value 
e res for trapped grains. It can also be treated as the maximum 
time of residence in the resonance for trapped grains. Finally, 
Tz.,, and T co u are the lifetimes of the resonant grains against 
collisions with background grains and other resonant grains, 
respectively. 

The actual lifetime of grains as determined by all sinks to- 
gether is given by 



1 1 

T Trf 



coll 



1 

Tcoll 



(7) 



and the maximum eccentricity achieved by the particles in the 
time T, as follows from Eq. @ for small e max , is 



&max — ^(-0 



q 



-BT. 



(8) 



p + q 

We return to Eq. (|6) and begin with its first term. The 
steady-state number of non-resonant grains in the annulus, no, 

is 



/;„ = /; ( l-^p)7> A >. 



(9) 



where the P-R drift time is given by (e.g. Burns et all 1 979 ) 

4fi max 



Tpr - 



B 



(10) 



Equation l|9} reflects the fact that the outer half of the annu- 
lus, the region between a(\ + e max ) and a, is always filled by 
the background grains, regardless of whether some of them get 
captured in the resonance upon arrival at a. In contrast, the 
number of grains in the inner half, between a and a(l — e max ), 
lies between zero (if p res ~\) and the number of grains in the 
outer half (if p res — 0). These two regions are shown in Fig. [2 
with dark and light grey, respectively. We note, however, that 
this division of the annulus into the outer and inner parts is done 
in terms of the semimajor axis and not in real space. Since grain 
orbits have non-zero eccentricities, no sharp boundary between 
the two halves exists in terms of distance, i.e. at r = a. The 
boundaries of the whole annulus will also be somewhat fuzzy. 

We introduce the normal optical depth To of the non- 
resonant population in the outer part of the annulus, which is 
not affected by the resonance. It is a product of the number of 
grains in the outer part and the cross section of a grain, divided 
by the area of the outer half of the annulus: 



to = 



n+T PR /2-ns 2 



Substituting Eq. dlOl and solving for n + results in 



n + =B{^)\o. 



(ID 



(12) 



Thus h + and thereby the first term in Eq. are expressed 
through an "observable" quantity, the background optical depth 

TO. 

The resonant pumping time T res in the second term of 
Eq. © comes from Eq. J4l >. For small e res , 



p + q. 



q B 



(13) 



Finally, the collisional lifetimes T®,, and T co u in Eq. (|6} 



should be estimated. General expressions are 



Vo 



' coll 



n a vo 



and 



Tcoll = 



V 
no v 



(14) 



where a = 4%s 2 is the collision cross section of the like-sized 
particles, vo the mean impact velocity between background and 
resonant grains, Vo the effective volume of their interaction 
specified below, and v and V are similar quantities for mutual 
collisions of resonant grains. 

Consider T® oll and the third term of Eq. (|6) fir st. The inter- 
action volume is given by (see, e.g., Appendix in lKrivov et al.l 
2006) 



(15) 



Vo = (2na) 2 h\ siny|. 



where h is the disk thickness (h = 2a sine for a disk with the 
semi-opening angle e) and y the angle between the orbits of the 
particles locked in the resonance and the orbits of non-resonant 
grains at their intersection point. With 



v*=(GM*(l-|3)/a) 



1/2 



(16) 



being the Keplerian circular velocity and neglecting corrections 
due to small eccentricities and inclinations, we have 



v w Vfc|sinY|. 



(17) 



Thus the ratio Vb/vo in the first Eq. (I14> is independent of ec- 
centricities and inclinations, as long as they are small, and that 
equation simplifies to 



' coll 



1 {2%a) 2 h 
n av k 



1 

no 



-Tcoll- 



(18) 



We also note that, in the monosize model considered, T^ is 
nearly independent of the particle size chosen. This is seen 
from Eqs. (|9) and <0, if 

we neglect a weak dependence of 

Vjt on the grain radius. 

To compute T co u for the fourth term of Eq. 0, we need 
the interaction volume V and the average impact velocity v for 
mutual collisions between the particles in the resonant clumps. 
Accurate calculations for this case are not easy, because the 
resonant population has a distribution of eccentricities from 
to e max and a highly non-uniform distribution of apsidal lines. 
There is no guarantee therefore that the particle-in-a-box-like 
estimates we employed for r% would yield correct results. A 
full solution of this problem, which may have a number of ap- 
plications to resonant systems, is beyond the scope of this pa- 
per. Our calculations (Queck et al., in prep.) show, however, 
that an analog of Eq. ( Tl8l . 



1 (2na) 2 h I. 

IcoU = — -Icolh 

n (JVjfc n 



(19) 



provides results that are within a factor of two from the "true" 
values. 

Thus all characteristic timescales, and all four terms on the 
right-hand side of Eq. (|6}, are now determined. A little techni- 
cal complication arises from the fact that T® in Eq. Jl 81 de- 
pends through no on e max (Eas.l9"UT0l which, in turn, depends 
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Fig. 2. Typical timescales as functions of the normal optical 
depth To of the non-resonant disk: Tpr (solid) T res (dashed), 
r% (dotted), T co n (dash-dotted). Thick, medium, and thin 
lines: p = 0.05, 0.10, and 0.20, respectively. These values cor- 
respond to s = 2.9, 1.5, and 0.7/jm. The collisional time T®,, 
is nearly the same for all three values (see text). 




[ - 

-7 -6 -5 -4 -3 

log (optical depth of the background disk) 

Fig. 3. Maximum orbital eccentricity of the clump particles as 
a function of the optical depth of the "regular" disk To- Dashed, 
solid, and dotted lines: e res = 0.3, 0.2, and 0.1, respectively. 
Thick, medium, and thin lines correspond to p = 0.05, 0.10, 
and 0.20. 



on T and thus back on T~L„ (Eqs.EHD- Actual calculations have 
therefore been done iteratively, with e res as the first approxima- 
tion for e max . 

For illustrative purposes, we choose the parameters suit- 
able for the e Erida ni disk: M* = . 8M^, the presumed 
planet at a n = 40 AU JLiouet alJl2000t lOzernov et aT]l200Ct 



lOuillen & Thorndikd l2002t iDeller & MaddisorJ l2005t) . The 
collisional lifetime Jl 8I > is proportional to the semi-opening 
angle of the disk £, which is unconstrained. By analogy with 
other debris disks seen edge-on (p Pic, AU Mic), we take e = 
0.1 = 6°. We choose 3:2-MMR, but the results are very simi- 
lar for other low-degree, low-order resonances. Figure [2]plots 
the timescales as functions of the optical depth To of the back- 
ground disk, showing that the collisional lifetime T® is shorter 
than the other timescales for disks with roughly To 10 -5 , 
which is the case for e Eri. 

Fi gure[3]shows how the maximum orbital eccentricity e max 
of resonant particles decreases with increasing optical depth. It 
also shows that, for denser disks, e max depends only weakly on 
the maximum value e res determined by the resonant dynamics. 
For To ~ 10~ 4 , the resulting value e max is close to 0.1, regard- 
less of e res . For such low eccentricities, the resonant popula- 
tion will look like a narrow ring without azimuthal structure, 
not like clumps. The dustier the disk, the narrower the result- 
ing ring (for planets in near-circular orbits), and the less pro- 
nounced the azimuthal structure. Therefore, prior to any cal- 
culation of the clumps' optical depth, we can conclude that in 
sufficiently dusty disks, scenario I fails to produce the clumps. 

We now solve the balance equation l|6). The collisional life- 
time of resonant particles against mutual collisions, Eq. dl9> . 
depends on n, which makes the whole equation Eq. (|6ji 



quadratic in n. At t 
cles tends to 



°°, the total number of resonant parti- 



TcoII 

2 



I4n+Pre 



oil 



r0 
1 coll 



Tre 



r0 

1 coll 



■(20) 



In the limit n <C no, meaning relatively faint clumps and a rela- 
tively bright "regular" disk, the fourth term of Eq. (|6) is much 
smaller than the third one, so that Eq. (|6j is linear in n. In this 
case, the above solution simplifies to 



n^Pres 



Tr-e} + {T? oll )- 



(21) 



It remains to estimate the optical depth of the resonant 
clumps. The clumps are not expected to be uniform for several 
reasons. First, even if the orbital eccentricities of all grains were 
equal, the second Kepler's law and the synodic motion of the 
particles relative to the planet would make the surface density 
of the resonant population non-uniform. Second, the resonant 
population consists of grains with a distribution of eccentric- 
ities, which is not uniform either: orbital dynamics forces the 
particle to spend most of the time in orbits with higher eccen- 
tricities up to e res (see. Eq.|4ji, whereas mutual collisions tend 
to eliminate these long-lived grains, shifting the distribution to- 
wards e <C e res . A superposition of these effects makes a spatial 
distribution of the optical depth in the clumps rather complex 
and dissimilar for different resonances and libration amplitudes 
A. 

To compute the normal optical depth T of the clumps, we 
use the following approximate method. We multiply n by the 
cross-section area of a particle, Ks 2 , and divide the result by 
an approximate area of the clumps S. Noting that the area of 
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the annulus is Sq = 4na 2 e max 
occupied by the clumps, S 



ASa 2 e 



we introduce the fractional area 
S/Sq. Then, 

(22) 



The fractional area S was determined with the aid of a straight- 
forward numerical simulation, by directly counting the points 
on the scatter plots. We considered the 3:2 resonance again, as- 
suming different values of e max between zero and e res = 0.2, 
as well as different values of A. For e max = 0.2, taking into 
account all spatial locations reached by the particles locked 
in the resonance leads to S from w 0.5 for A = 20° to ea 0.8 
for A = 90°. Defining, however, the regions with 30% of the 
peak brightness as a clump boundary, we arrive at values ss 0. 1 
for A = 20° to « 0.3 for A = 90°. Dependence on e max turned 
out to be rather weak: S grows slightly with decreasing emax- 
We therefore adopt a simple S = 0.2 in our numerical exam- 
ples, and stress that the actual x in the brightest portions of the 
clumps is several times larger than we predict. 

Figure|4](top) depicts the resulting optical depth of the res- 
onant blobs versus optical depth of the "regular" disk. The op- 
tical depth of clumps grows more slowly than that of the non- 
resonant population. A comparison between solid lines (Ea. l20t 
and dashed ones 

(Eq.ED of the same thickness shows that the 
linear approximation Mil provides very good accuracy. 

We now seek an approximate analytic scaling formula for 
the optical depth of the clumps, considering the stellar mass 
M*, the planet's orbital radius a p , and a number of other quan- 
tities as free parameters. The ratio of the normal optical depth 
of the clumps, x, and of the background disk, Xo, is 

_X_ _ n_ 1 —p res I '2 

To «o S 

From Eq. © and an obvious relation 

n = n + p res T, 

we obtain 
x T 



(23) 



(24) 



(25) 
T res and 
(26) 
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where Tpp is given by Eqs. dlOi and 0. 

For tenuous disks (such that T co ii 3> T res ), T 
e ma x(T) w e res - Therefore, 

T ^ *res "T re s 

— ~ Pres ■ - — Pres ' ~ 3 

Xo ST PR 4Se re s 

or, with the aid of Eq. < l 1 31 . 

X _ p + q e res 

— Pres ' ' , a • 

to q AS 

The same formula in a form that is convenient for numerical 
estimates is 



x 
to 



0125 (gj) 5 



S 



-l 



0.2 



p+q 
q 



(28) 



This equation shows that x/xo is typically smaller (but not 
much smaller) than unity. In some cases it can exceed unity, 
however. 
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Fig. 4. Normal optical depth of the resonant clumps x as a func- 
tion of the optical depth of the "regular" disk To. Top: full so- 
lution (I20> of the balance equation (solid lines) and solution 
(12 1 > of the linearized balance equation (dashed lines). Thick, 
medium, and thin lines correspond to p = 0.05, 0.10, and 0.20, 
respectively. Bottom: full solution J20l > of the balance equation 
for p = 0.1 (solid line) versus the tenuous disk approximation 
d27> (dotted) and denser disk approximation d29Kdashed). 

For denser disks (such that T co u <C T res ), the lifetime T in 
Eq. (1251 has to be evaluated iteratively, which complicates the 
calculations considerably. To keep the treatment at a reason- 
ably simple level, we therefore assume that n <C no, which is 
justified by Eq. d23l > and by the fact that x does not exceed Xo 
(27) for sufficiently dusty disks. In this case T sa T° = t C oll/ n 



coll 

We then substitute into Eq. \25\ expressions \\9\ for T co u, (|9} 
for no, d 1 Oft for Tpp, and so on "to the whole depth" to express 
them through the free parameters. Omitting straightforward, 
but lengthy algebra, the result is 

1 



Pres-~ § 



p7tsine 



c (l_fe) ( l_ P) 2/3 



IGM* 2 
— T o 



Krivov et al.: Clumps in debris disks 



7 



P + q 

, 1 



p + q 



1/3 



1 1/3 



(29) 



where the square-root expression is the planet's orbital velocity. 
Dropping factors in brackets that are close to unity, we obtain 
a somewhat cruder version of the same formula which is con- 
venient for numerical estimates: 

h / 1^21 ) ( J_\ 1 ( AS] ( L 

0.2/ V 0.1 / VO.lrad/ 



9x 10" 



,(Pres\ 
\0.5J 
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100AU 



-1/6 / T A 2 / 3 ( P </ 



2/3 



(30) 



The bottom panel in Fig.0]compares both analytic approxima- 
tions, Eqs. Mil and J29t . with each other and with the "exact" 
solution \2Q\ . 

Equations Mil and \29\ explicitly show the dependence of 
X on many parameters: the stellar mass M*, the planet's orbital 
radius a p , the resonance integers p and q, the typical p value of 
constituent particles, and others. For tenuous disks, x is nearly 

2/3 

proportional to Xrj. For denser disks, x °c Xq — a "slowdown 
effect seen before in Fig.|4] The clumps' optical depth is nearly 
independent of the planet's distance from the star a p for tenu- 
ous disks and slightly decreases with increasing a p for larger 
To. The dependence on the stellar mass M+ is absent for dilute 
disks and weak for denser ones: x slightly increases for more 
massive primaries. Note that the last statement is only true if 
we fix P and not s. If the grain radius s were treated as a free 
parameter, we would have to express (3 through s in Eq. i29\ . 
which would yield an additional in Eq. ( I29> . making the 
dependence on the stellar mass stronger. Next, an important de- 
pendence on the planet mass M p is implicit (through the capture 
probability p re s)\ we will discuss it later. The fractional area of 
the clumps, S, depends on the maximum orbital eccentricity 
of the resonant particles, and thus on other parameters such as 
To, but only weakly. For denser disks, x increases, albeit only 
slightly, with the semi-opening angle of the background disk. 

Figure [5] presents the "exact" solution, based on Eq. ( 1201 . 
in the form of contour plots of equal x, the average optical 
depth of the clumps. The figure shows that the phenomenon 
seen in Fig. |4] — "slowdown" of the growth of x at higher 



TO 



holds for all distances from the star a n . The effect is a 



bit more pronounced at larger a p . A comparison of the left and 
right panels demonstrates another effect: for given Xo and a p , 
disks around more massive/luminous stars can develop some- 
what denser clumps. All these effects are almost negligible, 
however. 

As for any analytic model, our model rests on a number 
of simplifying assumptions, which we now discuss. One of 
those is a constant value of the resonance trapping probabil- 
ity pres> which we have set to 0.5 in the numerical examples. 
In fact, it is through this quantity that the planet mass M p 
would enter the model. This probability also depends on a par- 
ticle's p, orbital semimajor axis a, and it varies from one res- 
onance to another. Non-zero inclinations can affect the trap- 
ping probability as well. To the best of our knowledge, ana- 
lytic or empirical formulas for p res as a function of M p are 
only available for some subspaces of these parameters, and 



only for t he planar problem (Beauge &Ferraz-Melio1ll994 
lLazzaro et alJI 19941) . It is clear, however, that resonance trap- 
ping is only possible within a certain planet mass range. If 
the plan et is too massive, the or bits of dust grains become 
chaotic (Kuchner& Holman 2003). For a solar-mass star, the 
maximum mass is typically several Jupiter masses. If, con- 
versely, the planet mass is too low, the particles go thro ugh 
the re sonance without getting trapped. Beauge & Ferraz-M ello 
ill 994 investigated the problem analytically and found, for in- 
stance, 0.04My (Mj = Jupiter mass) to be a minimum mass 
to trap gra ins with & = 01 in a 3:2 MMR. In a more gen- 
eral study, lLazzaro et alJ d 19941) found that for a planet with 
M p 2; O.OlMy and for micrometer-sized particles with eccen- 
tricities lower than a few percent, there are large regions of cap- 
ture in phase diagrams. These ana lytic estimates are supporte d 
by direct numerical integrations. iDeller & Maddisonl (|2005), 
for instance, found a high efficiency of capturing for a wide 
range of planet masses, from 0.01M/ to 3Mj, assuming a solar- 
mass central star. Within these limits p res shows rather a weak 
depende nce on M n , namely a moderate growth w ith M p /M* 
dBeauge & Ferraz-Mellol 1994 lLazzaro et al.ll994 . justifying 
our simple assumption p res = 0.5. 

Apart from a fixed p res , we used a fixed value e res of 0.2 
in the numerical examples. Contrary to p res , this quantity is 
not important at all for denser disks, because the maximum ec- 
centricity attained by the grains is limited by collisions, rather 
than by resonant dynamics; see Fig. [3] That is why e res does 
not enter Eqs. (I29i — d30i . For low optical depths, x is directly 
proportional to e res , see Eqs. d27l — d28i . As e res is known to lie 
approximately in the range [0.1,0.4], the uncertainty is only a 
factor of two even in this case. 

In the above treatment, we assumed destructive collisions 
to simply eliminate both colliders. This is supported by simple 
estimates. For two colliders of equal size s, the minimum rela- 
tive speed v needed to disrupt and disperse them, is determined 
by (see Eq. (5.2) of IKrivov et al.l2"ooi 



(31) 



where Q* D is the minimum specific energy for fragmentation 
and dispersal. A t dust sizes, Q* D ~ 10 7 to lO^rgg- 1 (e.g. 
Gri gorieva et alJ 12006). This gives v„„„ ~ 100 to 300ms _1 , 
much smaller than typical relative velocities (I17> . Under the 
same conditions, the largest fragment is an ord er of magni- 
tude s maller than either collider (see Sect. 5.3 of IKrivov et al.1 
2005), which allows us to ignore the generated fragments. We 
also ignore cratering collisions that may only moderately af- 
fect the mass budget, the timescales, and the optical depth of 
the clumps. 

Perhaps the most fundamental simplification that we make 
is to choose a single particle size. This leads, in particular, to 
overestimation of T®„ and T co u- The actual collisional lifetime 
will be shorter, because the grains that dominate the clumps 
will preferentially be destroyed by somewhat smaller, hence 
more abundant, particles. Collisions with interstellar grains 
may shorten the collisional lifetime even more. 

In the future, it would be interesting to develo p a more elab- 
orate kineti c model with n as a function of size s jRrivov et al.l 
2005, 2006). Such a model would be particularly interesting for 
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Fig. 5. Contour plots of the normal optical depth of the resonant clumps x, as a function of the optical depth of the "regular" disk 
To and the planet's distance from the star. Left and right: stars with 0.2 and 5.QM@. Contours are labeled with the values of log x. 
Assumed p = 0.10. 



dilute disks with Xo ^ 10 6 , expected to be resolved with future 
instruments (Herschel, ALMA, JWST, DARWIN, SAFIR, etc.) 



3. Scenario II 

In this section, we analyze the alternative scenario II, in which 
dust parent bodies are already locked in a certain resonance, 
so that the dust they produce resides in the same resonance, 
creating the clumps. 

This scenario is more difficult to quantify than the first one, 
for the following reasons. In scenario I, the dust production 
mechanisms were "hidden" in one quantity, the dust injection 
rate n + , which we estimated from an "observed" value, optical 
depth Xrj. This way essentially eliminates much of the uncer- 
tainty in modeling scenario I. In contrast, a physical descrip- 
tion of scenario II requires knowledge of various parameters of 
resonant planetesimal populations. These are not observable di- 
rectly and can only be accessed by extensive modeling, taking 
into account the wealth of processes that drive a forming plane- 
tary system (Wyatt 2003). Even if the parameters of a planetes- 
imal ensemble were known, another major source of uncertain- 
ties would be the necessity of modeling a collisional cascade 
from planetesimal-sized down to dust-sized objects, i.e. across 
some 30 orders of magnitude in mass. Therefore, we construct 
here the simplest possible model to explore the efficiency of 
scenario II, while being aware of its roughness. 

Consider a population of planetesimals with the total mass 
M, locked in the (p+q) : p- resonance. These planetesimals oc- 
casionally experience destructive collisions, in which smaller 
fragments are produced. Since the relative velocities of frag- 
ments are small compared to orbital velocities, most of them 
cannot leave the resonance and stay in the system. Gradual loss 
of material occurs through dust-sized fragments: below a cer- 
tain minimum radius s min , the radiation pressure will force the 
particles out of the resonance JWvatl20 06). Relative velocities, 



expected to be higher for smaller fragments, will also liberate 
a fraction of particles from the resonance. 

We numerically investigated the conditions under which the 
fragments remain captured in the resonance. Again, we took the 
e Eri system with = O.8M , M p = O.lMy, a p = 40 AU. We 
first chose the set of planetesimal orbits, all with the resonant 
(3:2) value of the semimajor axis, randomly chosen eccentrici- 
ties between and 0.1, inclinations between and 0.1 rad, and 
three angular elements between 0° and 360°, and selected those 
that turned out to be locked in resonance with libration ampli- 
tudes not larger than 30° (Fig.[6]top left). From the same set of 
planetesimals, we then released dust particles with different s 
(and P) and different velocity increments u relative to the par- 
ent bodies and numerically followed the motion of these ejecta 
to find out whether they will stay in the resonance as well. As 
a criterion for staying in the resonance, we simply require that 
the resonant argument <t> librates with a moderate amplitude A, 
not exceeding 30-50 degrees. We thus exclude shallow reso- 
nances that can easily be destroyed by external perturbations 
and, besides, lead to nearly axisymmetric ring-like configura- 
tions. 

The results are presented in Fig. [6] in the form of phase 
portraits of eccentricity and resonant angle, in polar coordi- 
nates. It is seen that particles larger than about 40/jm (or 
P < 0.003), if ejected with relative velocities u < 10ms , 
typically stay locked in the resonance. Thus p c ,, r «s 0.003 and 
Ucrit /vk ~ 0.003. The fact that p CT ,Y and u cr i t /vk are close to each 
other is not surprising. Both radiation pressure and an initial ve- 
locity "kick" do nothing else than cause the initial semimajor 
axis and eccentricity of the particle orbit to differ from those of 
a parent planetesimal. Within a factor of several, Aa/a ~ e ~ P 
in the first case and Aa/a ~ e ~ u/v^ in the second. Thus our 
finding means that a maximum relative change in (a, e) that still 
leaves the grains in the resonance is several tenths of a percent. 

Interestingly, a similar problem was investigated before 
in quite a different context. iRrivov & Banaszkiewiczl J200ll) 
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Fig. 6. Phase portraits of resonant planetesimals and the dust grains they release in mutual collisions. Top: planetesimals them- 
selves, middle: 40/jm-sized dust, bottom: 20/jm-sized dust. From left to right: the relative velocities of 0, 10, and 20ms~ 1 . 



studied the fate of tiny icy ejecta from the saturnian moon 
Hyperion, locked in the 4:3-MMR with Titan and found that 
about a half of 100/jm-sized grains remained locked in the res- 
onance. 

One of the two effects — the p-threshold — was re- 
cently considered in great detail bv lWvattl J2006h . For the 3:2- 
resonance, he found 



Per, 



: 0.034 



(Mp/MA 
\M*/M e J 



1/2 



(32) 



For M* = O.8M and M p = OAMj, this gives $ crit = 0.012, 
which is a slightly more "tolerant" a threshold than the one we 
observed in our simple simulation. 

As our numerical test was done for one pair (M t ,M p ), it 
still does not allow us to judge how the threshold of relative 
velocity scales with stellar mass and planetary mass. To find 
an approximate scaling rule, we therefore undertook a series 
of numerical integrations with 5 values of M p (0.03, 0.1, 0.3, 
1.0, and 3.0M/) and 6 values of M* (0.1, 0.2, 0.5, 1.0, 2.0, and 
5.0Mq), 5 x 6 = 30 runs in total. In each run, we integrated 
trajectories of 100 particles p = released with different ve- 
locities from planetesimals locked in the resonance with libra- 
tion amplitudes not larger than 30°. For each run separately, we 



determined the mean velocity threshold below which the par- 
ticles stay in the resonance, having libration amplitudes of not 
more than 30°, 40°, and 50° (Fig.0 symbols). We then fitted 
the results with a power law 



Ucrit/Vk = A 



( M p /Mj 
\M*/M e 



(33) 



The constants are A w 0.007 (±40%) and <B « 0.28 (±10%) for 
the "threshold" A of 30°. The fits are shown in Fig.Qwith lines. 
For the e Eri under the same assumptions as above, Eq. J33I 
gives u cr it = 14ms _1 for a planet with mass M p = O.lMj and 
26ms- 1 for M p = l.OMj. 

We note that the dust grains making the observed clumps 
have gone through a collisional cascade. The cascade starts 
with collisions of the largest planetesimals with mass m max (or 
radius s ma x), followed by collisions of progressively smaller 
bodies, until dust grains with mass m m ,„ (corresponds to s m ,„ 
or p CT ,Y explained above) are produced. We first consider a sin- 
gle collision and estimate the fraction of material which has 
u < U rrit- It can be approximated as jNakamura & Fuiiwaral 
11991b 



*P(< Ucrit) = 1 - ("0/ ' Ucrit) 1 (u c rit > «o) 



(34) 
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Fig. 7. Critical relative velocity of fragments u cr it/vk to stay 
in the resonance, as a function of the planet-to-star mass ratio. 
Symbols: results of individual runs; lines: power law fits (33). 
Different maximum A of particles were adopted as a criterion 
of staying in the resonance: A = 30° (triangles and solid line), 
A = 40° (squares and dashed line), A — 50° (stars and dotted 
line). 



where laboratory impact experiments suggest the lower cutoff 
mo ~5ms _1 (for dust sizes) and the exponent y ~ 2 (for con- 
solidated material, such as low-temperature ice). Of course, the 
power-law fit J34I should be treated with great caution. In real- 
ity, no sharp lower cutoff exists, so that up just places the lower 
bound on the speed range in which Eq. J34i is valid. A sin- 
gle power-law index y holds only in a limited velocity range. 
Both «o and y may vary significantly for different materials, 
morphology of the colliders, impact velocities, incident angles, 
targ et and proj ectile masses, and other parameters of collision 
(e.g. lFuiiwarall986l) . 

According to Eq. P4> . most of the fragments of a single 
collision have velocities slightly above mo. This is consistent 
with many impact experiments that repo rted relative velocities 
in the range 5-20 ms' 1 as typical ("e.g. lArakawa et al1ll995t 
iGiblinet al.l lT998). These typical velocities — and therefore 
the cutoff value up — d epend on the fragments' mass. Indeed, 
iNakamura & Fuiiwaral dl 99 ll) found a mass-velocity relation 
m °< m -1 / 6 to fit their experimental data. Since all impact ex- 
periments are only possible and have only been done for the 
lower end of the size distribution (objects less than one meter 
in size), we will assume that mo ~ 5 ms" 1 is valid for "the last" 
collisions in the cascade and that mo diminishes according to the 
m -1 / 6 law for "previous" collisions of larger objects. Needless 
to say, this prescription should also be treated with caution and 
is only suitable for crude estimates of the velocity-related ef- 
fects we are aiming at. 

A natural question is: what is the average number of colli- 
sions, k, needed to produce a dust particle with s„„„ = 30/jm out 
of parent planetesimals with size s max — 10 km? Simple analyt- 
ical calculations done in the Appendix yield k rj 8. This result 
depends, albeit rather weakly, on several parameters, such as 
the mass of the largest fragment in a single binary collision. 



We now estimate roughly the loss of particles due to u > 
u cr i t over the whole cascade. We simply replace the real cas- 
cade, in which each collision generates a mass distribution of 
fragments, with a sequence of k — 8 steps. In the first step, plan- 
etesimals of mass m max collide with each other and produce 
smaller fragments of equal mass. In the next step, the latter 
collide with each other, and so on, until the fragments of mass 
m m j n are produced. As follows from Eq. (1341 and the mass- 
velocity relation discussed above, the fraction of dust mass that 
is able to stay in the resonance can be estimated as 
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In nearly all situations, only the first term matters, which means 
that most of the material is lost to the "velocity effect" at the 
last collision. In what follows, we thus adopt 
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On the whole, the resonant population represents a nearly 
closed system, losing material only slowly through the smallest 
and fastest grains, which reduces the mass of the population on 
extremely long (typically Gyr) t imescales. This makes it pos- 
sible to apply Dohnanyi's ( 1969) theory of a closed collisional 
system. Provided that the critical impact energy for fragmen- 
tation and dispersal Q* D is independent of size over the size 
regime considered, he has shown that such a system evolves to 
a steady state, in which the mass distributio n is a simple power 
law: n (m] °c m^' a with a = 5/6 = 0.83. iDurda & DermotJ 
ill 9971) generalized this result by considering a more realistic, 
size-dependent Q* D . They have shown that a power-law depen- 
dence of Q* D on size preserves a power-law size distribution, 
but makes it steeper or flatter, depending on whether Q* D de- 
creases or increases with size. For Q* D °<- s~ 03 , which is a p lau- 
sible choice in the strength regime (see lKrivov et al . 2005, and 
references therein), Durda and Dermott found a = 0.87. 

The total cross section of the steady-state population is 
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(38) 



where p is the material density of the objects. For a = 5/6, 
Eq. J38l > takes a symmetric form 



(39) 
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Equations J38i and J39I show that the results depend on m m i„ 
and m mux only weakly. 

Like in the first scenario, we then divide E by the area of 
the clumps, S = 4nSa 2 e max , to obtain an estimate of the clumps' 
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Fig. 8. Normal optical depth of the resonant clumps X in sce- 
nario II, as a function of the planetesimal population mass 
M. Thick lines are for uq = 5ms" 1 and y = 2 and differ- 
ent a: Dohnanyi's a = 0.833 (dashed) and a more realistic 
a = 0.87 (solid). Thin lines are for a = 0.87 and different uq 
and y: uq = 1ms" 1 and y = 2 (solid), uq = 10ms" 1 and y = 2 
(dashed), uq = 5 ms" 1 and y = 1.5 (dotted). Cf. Fig.|4]for sce- 
nario I. 



normal optical depth. The results are presented in Fig.[8]for the 
e Eri case. We assumed again M* = O.8M , e = 0.1 =6°, and 
considered the 3:2 resonance (planetesimals with a = 56 AU) 
with a planet of mass M p = Q.lMj. We adopted s max — 10km, 
e,nax — 0.2, and S = 0.2. From ( l32l l. |3 6T ; r = 0.012 (or s,„,„ = 
12/jm). As in Fig.|4l we plot the normal optical depth of the 
clumps x, but now as a function of the planetesimal population 
mass M. Two size-distribution indexes were used: Dohnanyi's 
a = 0.833 and a more realistic a = 0.87. Our best-guess choice, 
a = 0.87 shows that the population with M = 0.01M e (M® 
is the Earth mass) would produce clumps with x ~ 10" 5 . The 
same figure demonstrates that the dependence on y is moderate, 
but that on the cutoff value uq is rather substantial. We specu- 
late that uncertainties in the fragment velocity distribution may 
introduce at least an order of magnitude, perhaps a larger, un- 
certainty in the resulting optical depth of the clumps. 

Similar to the first scenario, we now investigate the depen- 
dence of the clump formation efficiency on various parameters. 
We use Eq. J38l > with a = 0.87 and p = 2gcm" 3 and express 
m m i„ through P 6T ,- r and then through and M p by Eq. d32l >. By 
setting ^cascade t0 unity, we obtain a rough upper limit on the 
optical depth of the clumps: 




Each panel in Fig.|9]presents the results in form of the con- 
tour plots of equal x, the average optical depth of the clumps, 
on the total mass of the resonant planetesimals M and planet's 
orbital radius a p , showing that x grows with increasing M 
and decreasing a p . Different panels illustrate the dependence 
on a, M*, and M p . As already seen, the dependence on a is 
quite strong. The dependence on the stellar mass is strong, too. 
However, in contrast to the first scenario, stars with higher lu- 
minosity produce less pronounced clumps. Dependence on the 
planet mass M p is stronger than in the first scenario. More mas- 
sive planets create brighter clumps. 

Remember that all our estimates imply that the reso- 
nant population is in a steady state. Transient fluctuations 
in the popul ation can be produced by recent major colli- 
sional events ftVvatt & Dentll2002t iKenvon & Bromlevll2005t 
iGrigorieva et alJl2006l) . Usually, the dust debris cloud gener- 
ated by a "supercollision" of two large planetesimals is spread 
by the Keplerian shea r into a relatively ho mogeneous ring in 
only 10 4 to 10 5 years dKenvon & Bromlevl 120051) . However, in 
a resonant population, most of the fragments should stay in the 
resonance and contribute to the clump brightness — for the rea- 
sons discussed above. The same Eqs. J38i and J40l > apply, with 
M now the total mass of both colliders. Roughly, a collision 
of two lunar-sized objects would generate O.OlMfg of debris, 
enough to produce well-observable clumps. 

4. Conclusions and discussion 

In this paper, we address the possible origin of azimuthal sub- 
structure in the form of clumps or blobs, observed in some de- 
bris disks such as e Eridani. Two scenarios are considered. In 
the first case, the one discussed most in the literature, dust par- 
ticles are produced in planetesimal belts exterior to resonance 
locations, brought to exterior planetary resonances by the P-R 
effect, and their resonance trapping leads to the formation of 
clumps. In the second case, dust-producing planetesimals are 
already locked in a resonance, for instance, as a result of plan- 
etary migration in the past, and the dust they produce resides 
in the same resonance, creating the clumps. For both scenar- 
ios, we have constructed simple analytic models in order to ex- 
plore their applicability limits and efficiency for a wide range 
of stars, planets, disks, and planetesimal families. 

1. Scenario I. We conclude that the efficiency of the first 
scenario strongly depends on the normal geometrical optical 
depth Xo of the "regular" disk, from which dust is supplied 
to the clumps. The optical depth is known from observations. 
Besides, the brightness of the resulting clumps is nearly propor- 
tional to the probability of trapping "visible" dust grains in the 
resonance, p res . This parameter can be constrained relatively 
well, too, with the aid of numerical integrations for one or an- 
other particular system and a specified resonance. 

In disks with roughly Xo Ss 10" 4 -10" 5 , the first scenario 
works well, creating pronounced clumps. Their optical depth x 
is usually lower than, but can be comparable to or even exceed, 
the optical depth of the background disk. In other words, in di- 
lute disks the clumps can be almost as bright as, and sometimes 
brighter than, the underlying disk in the same region. This is 
the case in the solar system's interplanetary dust cloud and 
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Fig. 9. Contour plots of the normal optical depth of the resonant clumps x in scenario II, as a function of the planetesimal 
population mass M and the planet's distance from the star. Top left and top right: Dohnanyi's a = 0.833 and more realistic 
a = 0.87, both for the same star with O.8M and a planet with 0. 1 Jupiter mass. Middle left and middle right: stars with 0.5 and 
5.0Mq, both for a = 0.87 and a planet with 0. 1 Jupiter mass. Bottom left and bottom right: a planet with 0.03 Jupiter mass and 
1 Jupiter mass, both for a = 0.87 and M* = 0.8M (;) . Contours are labeled with the values of logx. Cf. Fig.|5]for scenario I. 
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the dust in the Kuiper belt with their To & 10~ 7 . Further, this 
will be the case for the debris disks at which future telescopes 
are aiming; examples are: the Herschel Space Observatory, the 
Atacama Large Millimeter Array (ALMA), the James Webb 
Space Telescope (JWST), the space infrared interferometer 
DARWIN, and the Spitzer's and Herschel's follow-on SAFIR. 
But it is not the case for debris disks of the other stars resolved 
so far. The disks currently observed and the clumps observed 
in them, all have optical depth > 10~ 5 . 

For higher optical depths To, X can still be high enough, 
but collisions are "killing" the clumps in the following way. 
They shorten the lifetimes of grains captured in the resonance, 
which prevents them from developing higher orbital eccentric- 
ities. Instead, the trapped particles remain in orbits with low 
eccentricities. The resonant population will look like a narrow 
ring without azimuthal structure, not like clumps. The dustier 
the disk, the narrower the resulting ring (for planets in near- 
circular orbits), and the less pronounced the azimuthal struc- 
ture. For dustier disks, we find that the optical depth of the 
clumps has a tendency to increase when the responsible planet 
orbits closer to the star and has a higher mass. It also tends 
to grow with the luminosity of the central star. Increasing the 
opening angle of the dust disk has the same effect. However, 
all these dependencies are quite weak. 

In the case of e Eridani, which falls into the category of 
"dustier" disks, the first scenario does not seem to be consistent 
with the presence of several longitudinally confined "blobs". 

2. Scenario II. The efficiency of the second mechanism is 
obviously proportional to the total mass of the resonant popula- 
tion of planetesimals, but is strongly affected by other parame- 
ters, too. The brightness of the clumps produced by the second 
scenario increases with the decreasing luminosity of the star, 
increasing planetary mass, and decreasing orbital radius of the 
planet. All these dependencies are more pronounced than in the 
first scenario. 

In the case of e Eridani, a family with a total mass of <~ 0.01 
to 0.1 Earth masses could well be sufficient to account for the 
observed optical depth of the clumps, roug hly T - 10~ 5 . Thus 
this mechanism appears to be more relevant for systems like 
e Eri. 

It must be noted that the efficiency of scenario II is quite 
sensitive to the poorly known properties of the collisional 
grinding process. One is the critical energy for fragmentation 
and dispersal Q* D over the whole range of masses — from 
planetesimals down to dust — which strongly affects the equi- 
librium slope of the size distribution, and therefore the total 
cross section ares of the resonant clumps. Even more impor- 
tant is the distribution of the ejecta velocities in disruptive col- 
lisions, which determines what fraction of collisional debris is 
fast enough to leave the resonant population. For these reasons, 
models of the second scenario are quantitatively more uncer- 
tain than those of the first one. 

Models presented here are the simplest possible ones we 
could construct. Accordingly, they involve a wealth of simpli- 
fying assumptions and therefore cannot substitute for develop- 
ing more elaborate models. Yet we hope that our results can 
give useful guidelines and serve as a starting point for more 
detailed studies of clumps in debris disks. 
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Appendix A: Number of collisions over the 
collisional cascade 

Consider the collisional cascade that produces small grains of mass 
m from parent bodies of mass m max . We wish to estimate the mean 
number of collisional events that the grains of mass m must have had 
in their collisional history. 

Each individual fragmenting collision of this cascade can be de- 
scribed by a distribution g(m,m x ) of fragments m from the largest one 
m x down to infinitely small ones. Let the normalization condition be 
given by 



/ 



g(m,m x )dm = 1. 



Now we introduce a relation between m x and m max given by /j = 
m x /m max . As each catastrophic collision in a cascade reduces the mass 
of the largest possible fragment by the factor fj, the maximum number 
of collisions that a fragment of size m might have undergone since the 
parent body m max was disrupted is given by 

_ ln(m/m max ) 
hijU 

If we assume every collision to be described by the same fj and the 
same power-law dependence for 

ti-2 

i \ i \ m x 
g(m,m x ) = (2-n)- ?Fr , 

but with a different absolute value of m x , the resulting distribution 
after k + 1 collisions is given by 



mi/i 



m 2 fi 



P k {m,m max ) = j g(mi,m max v) J g(m 2 ,mnj) J g(m 3 ,m 2 lj) ■ 



j g(m k ,m k -ifi) g{m,m k n) 



dm k ■ ■ ■ dm 3 dm2dmi , 



(A.1) 



where P k provides the fraction of mass that falls in the range [m,m + 
dm] after k + 1 collisions. If we change to logarithmic scaling of the 
integration variables, i.e. dm\ = mid(ln^J-) and so on, and explicitly 
insert g, we obtain 



Pk{m,m max ) = (2-r]) k+1 



tl-1 / Tl — 1 

ln/i-* 1 In/i 1 -* 2 



/ 



(myi) 



ti-2 



/ 



m\ 1 J m) 1 m 1 ! 1 

2-t 3 In//- 2 k 



x m k ■ ■ ■ m^m2mi 

d(ln — ) ■ • ■ d(ln — )d(ln — )d(ln — 
m m m m 



(A.2) 
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Luckily, all integration variables disappear and the integrand is con- 
stant: 



0.07 



Pk(m,m max ) = (2-T|) 
In 



/ k+l \ 'H -2 

fc+1 / P m max 



m \ m 

/ / / - / 



tofr k \n^- k ln/j 2 -* 



Ah m k\ Ah m T>\Ah m 2\Ah 

d(ln — ) • ■ • d(ln — )d(ln — )d(ln — ) 
m m m m 

Now the variables are shifted starting with m\ — > f/m\ : 
n*+1 1 ( 



(A.3) 



P k {m,m max ) = (2-T))* 



7)7 



'max 



m 
In =2 



/ / /-/ 



lnp 1 



ln/j 2 



ln/i- 2 



j/i -K\ All m ?>\Ah m 2\Ah m l 

d(ln — ) • ■ • d(ln — )d(ln — )d(ln — 
m m m m 



(A.4) 



Substituting everything from m\ to m k in the same manner, we obtain 

n -2 



P k (m,m„ 



(2 -11) 



1 fn k+x m„ 



I I I- 1 





Ah Ah m ^\Ah m ^\Ah '"1 

d(ln — ) • ■ • d(ln — )d(ln — )d(ln — 
m m m m 

The integral equates the volume spanned by < x k < ■ ■ 
x\ <xq = In - — ■^ !m , which is given by 



(A.5) 

< X3 < X2 < 



-ro=ln ' 

/ 





dxic ■ ■ -dx3dx2dxi 





k 



k\ 

The final result is 

Pk(m,m max ) 



m \ m 



(A.6) 



n-2 



k\ 



(A.7) 



Pk(m,m max ) 
2-T] 



1-2 ^ In 



2-n 



(A.8) 



ra \ m J k\ 
For & = 0, i.e. no intermediate collision, we are back at 

P (m,m max ) = g(m,^m max ). (A.9) 

To obtain the distribution of the total mass of fragments per mass 
decade, we have to evaluate 

P = mln(10)/>. (A. 10) 

Figure lA~T1 shows examples for the distribution. The planetesimal ra- 
dius s max = 10km and the grain radius s = 100 to 10/jm correspond 
to the mass ratio log l0 (m max /m) ~ 24 to 27, or x from —3 to 0. In this 
case, the distributions peak at k = 7 to 10. Therefore, we adopt k = 8 
for the bulk of the calculations. 




6 8 10 

Number of collisions 

Fig. A.l. Probability distribution of particles with size m orig- 
inating from parent bodies of mass m max = 10 27+A m over the 
number of collisions they underwent, x has values —6, —3, 0, 
3, 6 (peaks from left to right). r| = 11/6. The peak position 
is roughly proportional to log 10 (m muv /m)/3.5. (The lines just 
connect adjacent points for illustrative purpose, although the 
distribution is discrete in k. The distribution looks similar to a 
Poisson distribution but is different.) Top: m x /m max — jj = 0.1, 
bottom: /j = 0.5. 
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